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Abstract 

The MSSM with exphcit CP violation is studied through the di-photon decay channel of the 
lightest neutral Higgs boson. Through the leading one- loop order Hi 77 is affected by a large 
number of Higgs-sparticle couplings, which could be complex. Our preliminary scan over the 
Supersymmetric parameter space shows that more than 50% average deviations are possible, in 
either direction, in the corresponding branching ratio, with respect to the case of the CP-conserving 
MSSM. In particular, our analysis shows that in the presence of a light stop (with mass ^ 200 
GeV) a CP-violating phase (j)^ ~ 90° can render the Hi — > 77 branching ratio more than 10 times 
larger, for suitable combinations of the other MSSM parameters. 
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The mechanism of Electro- Weak Symmetry Breaking (EWSB) is elusive even after the very 
successful LEP era, although precision measurements hint at a light Higgs particle. It is 
expected that the soon to be operational Large Hadron Collider (LHC) will be able to 
make definite statements about the Higgs mechanism. At the same time there are various 
reasons to think that the Standard Model (SM) is only an effective theory valid up to TeV 
range, and some richer structure is needed to explain particle dynamics (much) beyond 
such energy scale. Supersymmetry (SUSY), being the most favoured of all the new physics 
scenarios proposed so far, is going to be searched for in all possible ways at the LHC. Within 
the Minimal Supersymmetric Standard Model (MSSM) the scalar potential conserves CP at 
tree level [l]. The reason is that SUSY imposes an additional (holomorphic) symmetry on 
the Higgs sector of a general two-Higgs doublet model, that entails flavour conservation in 
tree-level neutral currents and absence of CP-violating scalar-pseudoscalar mixings in the 
Born approximation. Beyond the latter, recent studies have shown that CP invariance of 
the Higgs potential may in principle be broken by radiative corrections [2|], as the Vacuum 
Expectation Values ( VEVs) of the two Higgs doublets can develop a relative phase . This 
type of CP violation is generally referred to as spontaneous CP violation and it requires a 
light Higgs state as a result of the Georg-Pa-s theorem g but the poss.bihty of the latter 
has now been essentially ruled out by experiment [5[. 

CP violation can also be explicitly induced in the MSSM, in much the same way as it 
is done in the SM, by complex Yukawa couplings of the Higgs bosons to (s)quarks. There 
are several new parameters in the Supersymmetric theory, that are absent in the SM, which 
could well be complex and thus possess CP-violating phases. Such parameters include: (i) 
the higgsino mass term /i, (ii) the soft SUSY-breaking gaugino masses Ma (a = 1,2,3), 
(iii) the soft bilinear term Bfi and (iv) the soft trilinear Yukawa couplings Af of the Higgs 
particles to scalar fermions of flavour /. 

Each of these parameters can have independent phases. After applying universality con- 
ditions at a unification scale Mx the gaugino masses have a common phase and the trilinear 
couplings are all equal with another common phase. As argued by [6], one may deviate 
from exact universality and consider A f to he diagonal in flavour space with vanishing first 
and second generation couplings to avoid problems with the electron, muon and neutron 
Electric Dipole Moments (EDMs). This leaves four independent phases, those of fi,Bfi,Ma 
and Af. However, the two U{1) symmetries of the conformal-invariant part of the MSSM 
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may be employed to re-phase one of the Higgs doublet fields and the gaugino fields such 
that Ma and Bjj are real 0, 0]. We will work within this setup with two independent 
physical phases, which we take to be arg{fi) = 0^ and arg{Af) = (p^^.. As intimated, the 
CP- violating phases associated with the sfermions of the first and, to a lesser extent, sec- 
ond generations are severely constrained by bounds on the EDMs of the electron, neutron 
and muon. However, there have been several suggestions to evade these constraints 



without suppressing the CP-violating phases. One possibility is to arrange for partial can- 
cellations among various contributions to the EDMs ll| . Another option is to make the first 
two generations of scalar fermions rather heavy, of order a few TeV, so that the one-loop 
EDM constraints are automatically evaded. As a matter of fact, one can consider so-called 



effective SUSY models 10| where decoupling of the first and second generation sfermions are 
invoked to solve the SUSY Flavour Changing Neutral Current (FCNC) and CP problems 
without spoiling the naturalness condition. We adopted the latter version of a CP-violating 
MSSM for our analysis, along with Af = for the first two generation sfermions. 

The CP- violating phases (p^ and (pAf could in principle be measured directly in the pro; 



duction cross sections and decay widths of (s)particles in high energy colliders j?], 12|] - 17 1 
or indirectly via their radiative effects on the Higgs sector p, [l^ . In this letter we will look 
at Hi — > 77 which involve the (leading) direct effects of CP violation through couplings of 
Hi to sparticles in the loops (see Fig. [1]) as well as the (subleading) indirect effect through 
the scalar-pseudoscalar mixing yielding a CP-mixed state. Hi. The origin of this CP-mixing 
is the following. In the Higgs sector, the CP-violating phases mentioned above introduce 
non-vanishing off-diagonal mixing terms in the neutral Higgs mass matrix, which in the 
weak basis (0i,02,o), when 
schematically be written as 



weak basis (0i,02;O), where 0i2 are the CP-even states and a is the CP-odd state, may 

B,QQ,Q 

, ( Ml Mlp\ 
[Mis Ml ) 

Here, is a 2 x 2 matrix describing the transition between the CP-even states, Ai^, 
gives the mass of the CP-odd state whilst A^p^ = (7Vl|p)^ (a 1 x 2 matrix) describes the 
mixing between the CP-even and CP-odd states. The mixing matrix elements are typically 
proportional to 

MIp oc Im{^Af) (2) 



and are dominated by loops involving the top squarks and could be of order M|. As a result, 
the neutral Higgs bosons of the MSSM no longer carry any definite CP-parities. Rotation 
from the EW states to the mass eigenvalues, 

is now carried out by a 3 x 3 real orthogonal matrix O, such that 



O^A^^O = diag(M^^, M^^, M^J 



(3) 



with Mhj^ < < Mh-^- As a consequence, it is now appropriate to parameterise the Higgs 
sector of the CP- violating MSSM in terms of the mass of the charged Higgs boson, Mh±, 
as the latter remains basically unaffected. (For a detailed formulation of the MSSM Higgs 



sector with explicit CP violation, see Refs. 



3, 



18|.) 



In order to study the effects of the CP-violating phases we focus here on the di-photon 
decay mode of the lightest neutral Higgs boson. Hi. The reason is twofold. Firstly, the 
di-photon decay mode is the most promising channel for the discovery of a light neutral 
Higgs state - of mass between, say, 80 — 130 GeV, at the LHC 
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2l| . Secondly, the 



dominant CP- violating terms dependent on fi and Af (hereafter, / = b,t,T) enter the 
perturbative calculation of the di-photon decay width with a coupling strength that is of 
the same order as that of the CP-conserving ones (of 0{a^)). Furthermore, on the technical 
side, thanks to the narrow width of such a light Higgs state (of 10 MeV at the most), the 
entire gg/qq ^ Hi ^ 77 process can be factorised into three parts: the production process, 
the Higgs propagator and the decay channel. Effects of CP violation can show up in this 
process through the aforementioned couplings in the production, through a possible mixing 
of Higgs states at one-loop and above in the propagator and through the same couplings 
in the decay. CP violation entering the production of a Higgs state in gluon-gluon fusion 
process at hadron colliders was studied first by I22I, choosing a parameter space region 



which is not sensitive to the CP-mixing of the Higgs states, and later by 
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23|, including 



the presence of CP-mixing of the Higgs states. Effects of CP-mixing in the propagator are 



discussed separately but in great detail in pi3|]. A thorough study o: 



the other MSSM Higgs 



decay channels in presence of CP-violation can be found in [l2]-|18|],[22|-[25(]. We postpone 
the full analyses of gg/qq Hi 'j'j including also CP-violating effects in production and 
propagation to a later work. 
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FIG. 1: Diagrams for Higgs decay into 77 pairs in the CP-violating MSSM: f = t, b; f = 

In the MSSM a Higgs state decays into two photons through loops of fermions, sfermions, 
EW gauge bosons, charged Higgses as well as charginos (see Fig. [1]). (Expressions for the 
amplitudes of Hi 77 along with relevant couplings are available in [26| and references 
therein.) This decay mode of if 1 (also H2 and H3) along with its production through gluon- 



gluon fusion is discussed by Ref . 23| . However, the analysis therein was limited to parameter 
space regions where CP-violating effects are only due to the changed SM particle (especially 
W^) couplings to the Hi. Effects of sparticles (made suitably heavy in Ref. (23|) in the 
loops were negligible. We examine here the complementary region of MSSM parameter 
space (albeit limitedly to the Higgs decay), wherein sparticle masses are taken light, so that 
they contribute substantially in the loops. In addition, we will show that, in the presence 
of non-trivial CP-violating phases, there are regions of MSSM parameter space where the 
couplings of the decaying Higgs bosons to all sparticles in the loops are strongly modified with 
respect to the CP-conserving MSSM, thereby inducing dramatic changes on the Hi '-fj 
width and Branching Ratio (BR). 

n 

To prove this, we have used the publicly available FORTRAN code CPSuperH [29] version 
2 for our analysis. CPSuperH calculates the mass spectrum and decay widths of the neu- 
tral and charged Higgs bosons in the most general MSSM including explicitly CP-violating 
phases. In addition, it computes all the couplings of the neutral Higgs bosons Hi 2,3 and the 
charged Higgs boson if ^ to ordinary and Supersymmetric matter. The program is based on 
the results obtained in Refs. [3]-[l5| and the most recent renormalisation group improved 
effective-potential approach, which includes dominant higher-order logarithmic and thresh- 
old corrections, 6-quark Yukawa-coupling resummation effects and Higgs boson pole-mass 



shifts 
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27|. 
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The free non-SM parameters of the model now include: phase of /i (0^), charged 
Higgs mass {Mh±), soft gaugino masses {Ma), soft sfermion masses of the third generation 
{M(Q.^ Us D3 is Es))^ (unified) soft trilinear couplings of the third generation (|A/|), phase 
of the trilinear coupling {(pAf)- In om scan we have chosen the following very extensive 
parameter ranges: 

tan/3 : 1 - 60, : 100 - 2000 GeV, (f)^, : 0° - 180°, Mh± : 100 - 400 GeV, 
M2 : 100 - 500 GeV, M^g^ ^^^^^^^^^^^^^ : 100 - 2000 GeV, \Af \ : 100 - 2000 GeV. 



We aimed at searching regions in the Supersymmetric parameter space where the variation 
in the BR(iyi 77) due to the CP-violating phases is maximised compared to the CP- 
conserving case. As stated above, the CP- violating effects are proportional to Xm(yuAj), 
and so we opted to fix (pAf to 0° and varied only (p^. Besides, Mi and M3 were kept fixed 
as their variation is of no significance here. Finally, no unification of the soft sfermion 
masses was assumed, though they were taken in the same range. For this analysis, threshold 
corrections induced by the exchange of gluinos and charginos in the Higgs-quark-antiquark 



vertices 



28 



291] were not included. 

We scanned the above parameter space for 100,000 randomly selected points within the 
ranges specified above and for each of these we have taken values of 0^ increasing from 
0° to 180° in steps of 20°. Notice that 0^ = 0°(180)° corresponds to the CP-conserving 
MSSM point with /i = while any other non-trivial 0^ shows the effect of CP 



violation. The following experimental constraints from LEP2 and Tevatron j30l . l3l| were 
imposed during the scan: 

m^± > 104 GeV (LEP2), 
m? > 100 GeV for / = [, i>, k (LEP2), 



> 300 GeV (Tevatron). 



Parameter space points violating these constraints were discarded. We only considered 
points with lightest Higgs mass (Mhi) between 90 and 130 GeV, the range in which the 
Hi — > 77 decay is relevant. 

For each point in the scans that survives the various constraints discussed above we asked 
CPSuperH to print out the mass and 77 BR of the lightest Higgs Hi. In order to have an 



idea of the overall trend followed by the BR for different phases, we first looked at the 
average behaviour at specific M//^ values. To do this we divided the mass range into bins of 
size 4 GeV. To find the average sensitivity within each mass bin, we defined the percentage 
deviation 

^ J ^ 



where the summation is over the number of random points {n) within a particular bin i. We 
denote the BR in the CP-conserving case by BRq (specifically, the latter corresponds to the 
case 0^ = 0°, however, without any loss or gain of information, we could alternatively have 
used the limit 0^ = 180°) and that with a non-vanishing 0^ (different from 180°) by BR^. 
This average percentage deviation is plotted in Fig. [2] (left) for the different values of 0^ 
taken in each bin. There is an enhancement in the BR of about 20% for Mh^ larger than 110 
GeV for moderate values of (p^ ~ 100°, while there is a suppression of about 40 — 50% for 
Mh^ around 90 - 98 GeV. This change-over from enhancement to suppression for lower Mh^ 
values shows the diminishing role of sparticles in the loop as the mass difference between 
2m j and M^^ increases and the effect of a non-zero 0^ is effectively more and more through 



a changed HiWW coupling. Such suppression is in agreement with the results of [23|]. In 
the mass region of 100 - 110 GeV the effect is apparently very small. However, this is an 
artifact of the binned averaging, where points with enhanced and suppressed BRs falling in 
the same mass bin cancel each other. This cancellation is nullified by taking the absolute 
value of the difference in the numerator of Eq. (jlj), i.e., before averaging. The result is 
plotted in Fig. [2] (right). More than 50% deviation is seen for 0^ = 100° for M^^ around 
104 GeV. Now, it should be noted that these figures represent only the average behaviour. 
It is therefore possible to find regions of parameter space where the differences are larger (or 
smaller, for that matter). We did indeed find points with difference in the BR larger than 
50 times in our scan in either direction. 

A subtlety should be noted in this context though, concerning the derived MSSM masses 
that also depend on 0^ and enter the decay Hi 77 (M/^-^, mi^-, M^±). In fact, all the 
latter change when going from the CP-conserving case to the CP-violating one. The most 
crucial one in this respect is Mh^- However, we have verified that for the same parameter 
point (apart from a different 0^) the latter always changes less than 2 GeV between the two 
MSSM configurations. Hence, our 4 GeV wide bins do capture percentage corrections to the 
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FIG. 2: Scan result showing binned average values of percentage differences in the BR(ifi 77) 
between the CP-violating and CP-conserving case (left) as well as the absolute value of it (right) 
- see Eq. ([4]) - for various choices of 0^. 



BR consistently between the two MSSMs as a function of the Hghtest Higgs boson mass. 
(Rare borderline cases are also correctly assigned to the right bin.) Besides, 2 GeV is roughly 
the di-photon mass resolution in ATLAS [2C| (while in CMS it is somewhat better 2l|). In 



short, we imagine an experimental situation in which a Higgs resonance is extracted in 77X 
samples with the above mass resolution at a time when the other SUSY masses and mixing 
(including tan/3) entering the loops of the Hi 77 mode have already been measured in 
real sparticle production with a resolution that does not allow one to distinguish between a 
CP-conserving and a CP-violating MSSM scenario. Under these conditions, for large enough 
differences of BRs, a simple measurement of the normalisation of the 77 resonance (after 
background subtraction) may suffice to distinguish between the two envisaged CP scenarios. 



To illustrate the validity of this argument, we have selected some specific points from 
the parameter scan to study the deviations in the BR. While doing this we also wanted to 
understand the contribution of the different components inducing the CP- violating effect. 
We specifically looked for points with large trilinear coupling and large fi values, for the CP- 
mixing of the Higgs states is proportional to their product (as mentioned in the beginning 
of this letter). We then considered different cases choosing the soft-mass values such that 
only one of the sparticles in the loop is light, while all others are heavy (the mass of the 
charged Higges boson is varied between 100 GeV and 400 GeV, but for the interesting region 
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FIG. 3: BR(i?i 77) plotted against M^i- Parameters used are: tan/3 = 20, Mi = 
100 GeV, M2 = M3 = 1 TeV, ^((33,03,13,^3) = 1 TeV, = 1 TeV, \Af\ = 1.5 TeV. 

Left figure has = 1 TeV while right one has M^^^ = 250 GeV (the latter giving a rather light 
stop, = 200 GeV). 



> 115 GeV it is heavier than 300 GeV), expecting to see the effect due only to the 
exchange of this light sparticle, alongside the one due to standard matter, t, b and W^. 
We also considered the situation when only SM particles are effective, with all the relevant 
sparticles heavy. We plot the BR against M^i for 0^ = 0°, 90°, 180° for two cases, (i) with 
all sparticles heavy and (ii) with a light stop of around 200 GeV, in Fig. [31 In case (i) the 
effect is almost entirely due to the CP-mixing of the Higgs states, entering the BR through 
the deviation in the couplings of the SM particles with the Hi. In case (ii) we also have, 
in addition to the above, the effect of a light stop, through its CP-violating coupling with 
the Hi as well as the sensitivity of the stop mass to the CP-violating phases. It is instead 
found that the effects of light sbottoms, staus and charginos are negligible, so that the BRs 
in these cases - keeping all other SUSY parameters to be the same - are similar to case (i). 
This is indeed expected, considering the smaller Yukawa couplings for the corresponding SM 
particles (with respect to the top) and - in particular - the stringent experimental limit on 
m^^ from Tevatron. The effect in the case of only heavy sparticles in the loop is exclusively 
due to the modification of the SM fermion and gauge boson couplings of the CP-mixed Higgs 
state. 

As stated above, the deviations in the Mh^ is within 2 GeV for both the cases for the 
entire range of 0^, as shown in Fig. H] (left). The larger sensitivity of the BR to 0^ in the case 
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FIG. 4: Mff-^ (left) and m^^ (right) plotted against (p^^. In the case of Mfji the solid line corresponds 
to the case with Mjj^ = 1 TeV while the dotted line corresponds to Mjj_^ = 250 GeV. Other SUSY 
parameters are the ones given in Fig. [3l 

with light stop can partly be explained through the sensitivity of m^^ to 0^. The variation 
of m^-^ with 0^ is plotted in Fig. |4] (right). Notice that the latter is of the same order as 



the expected experimental resolution 20|, l2l||, so that it may not be possible to confirm 



CP-violating effects directly in the stop sector. 

In summary, while a full study incorporating the production processes and detector de- 
pendent aspects is needed to have a clear quantitative picture, our prehminary analyses in- 
dicate that the di-photon channel of the lightest Higgs boson may enable one to distinguish 
the CP-violating MSSM from the CP-conserving one, so long that some SUSY parameters 
are measured elsewhere. This is not phenomenologically unconceivable, as the Hi 77 
detection mode requires a very high luminosity, unlike the discovery of those sparticles (and 
the measurement of their masses and couplings) that impinge on the Higgs process studied 
here. A complete analysis will eventually require to fold the decay process with propagator 
effects and the appropriate production mode (gluon-gluon fusion and Higgs-strahlung in this 
case), where similar CP- violating effects may enter. 
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